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Abstract. Multiple orthogonal polynomials are traditionally studied because 
of their connections to number theory and approximation theory. In recent 
years they were found to be connected to certain models in random matrix 
theory. In this paper we introduce the notion of a multiple orthogonal polyno- 
mial ensemble (MOP ensemble) and derive some of their basic properties. It is 
shown that Angelesco and Nikishin systems give rise to MOP ensembles and 
that the equilibrium problems that are associated with these systems have a 
natural interpretation in the context of MOP ensembles. 



1. Introduction 

Multiple Orthogonal Polynomials (MOPs) were introduced and studied for 
problems in analytic number theory (irrationality and transcendence proofs). Later 
they appeared in approximation theory, most notably in the theory of Hermite-Pade 
approximation and in this context they are also called Hermite-Pade polynomials 
[aiilllOllMllMlliQlllllliailllleO]. MOPs were also studied from the point of 
view of new special functions [fil [HI [HI |49l [Ml IZ] . See the books [Ml [5^ and 
the survey papers [3, 8, 165|, 166] for these aspects of MOPs. Further developments 
in these directions are reported in e.g. [71 [TUl [TBI [2TI [23l [24l [371 l38l [50l l54j . 

Recently MOPs also appeared in a natural way in probability theory and 
mathematical physics in certain models coming from random matrix theory and 
non- intersecting paths. The connection was first observed in [14] where MOPs 
were used in in a random matrix model with external source. In the Gauss- 
ian case, the external source model has an equivalent interpretation in terms of 
non-intersecting Brownian motions. The external source model was further an- 
alyzed with the use of multiple Hermite and multiple Laguerre polynomials in 
[5l[9l[l5l[lll[3ll[43l[6ll[53l[55l [Ml [69], see also [S [SI [IS [30l [M]. A related 
non-intersecting path model was studied in |46j using MOPs for modified Bessel 
weights that were introduced earlier in [22] . The biorthogonal polynomials aris- 
ing in the two matrix model were identified as MOPs in |47) . For a special case 
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they were asymptotically analyzed in |36|, I57j . The Cauchy two matrix model and 
their associated Cauchy biorthogonal polynomials have a number of similar features 
[12l [T3] . MOPs were generelized to MOPs of mixed type in [H HH [271128] , 

Asymptotic results were mainly obtained from an analysis of the Riemann- 
Hilbert problem for MOPs, formulated by Van Assche et al. j68j as an extension of 
the Riemann-Hilbert problem for orthogonal polynomials [39 . The application of 
the Deift/Zhou steepest descent analysis [29 to the Riemann-Hilbert problem for 
MOPs presents several interesting new features that however we will not discuss 
here. 

It is the aim of this paper to give an introductory account of MOPs from the 
point of view of determinantal point processes. After discussing the definition and 
some of the basic properties of MOPs we discuss a multiple integral representatons 
for the type II MOPs, which is essentially taken from |14j . Under a suitable constant 
sign condition the formula can be interpreted as the expectation value of the random 
polynomial 11^^]^ (2; — Xj) with roots xi, . . . , a;„ from a determinantal point process 
(called a MOP ensemble) on the real line. 

The constant sign condition holds in particular for Angelesco and Nikishin 
systems. For both of these systems we show that the joint p.d.f. of the associated 
MOP ensemble takes on a particular nice form. In the large n limit it allows for 
a natural probabilistic interpretation of the vector equilibrium problems that are 
associated with Angelesco and Nikishin systems. 



2.1. Definitions. Given weight functions wi, . . . , Wp on M and a multi-index 
n — {ni,...,np) G N^, the type II MOP is a monic polynomial Pft of degree 
\n\ = ni + ■ ■ ■ + Hp such that 



The conditions (|2.ip give a system of n linear equations for the n free coefficients 
of the polynomial Pfi (recall that is monic) . If the system has a unique solution 
we say that the multi- index n is normal (with respect to the weights wi, . . . , Wp). 

In this paper we mainly deal with the type 11 MOP, but at times it is useful 

(i) 

to consider the dual notion of type I MOPs as well. These are polynomials , 
j — 1, . . . ,p, of degrees deg A^-* — nj — I, such that the linear form 



2. Multiple orthogonal polynomials 



(2.1) 




k = 




Throughout we will write 



n — \n\ = ni + ■ ■ ■ + n,^ 



p 



(2.2) 




satisfies 



(2.3) 




k = 



0,l,...,n-2. 



If we supplement this with the normalizing condition 



(2.4) 
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then again we have a system of n = |n| hnear equations for the in total n coefficients 
of the polynomials A'^^\ j — 1, . . . ,p. 



2.2. Determinantal expressions. Let 



„(i) 



x'^Wj {x)dx 



denote the fcth moment of the weight Wj , and let 

rr(i) _ (M \ 



k—0,...,m,l—0,...,n 



be the (m+l) x (n+1) Hankel matrix with the moments of Wj. The conditions (|2.3p 
and (|2.4p give rise to a linear system whose matrix has the block Hankel structure 



(2.5) 



Ma = 



rr(l) rr(2) 

?i— l.ni — 1 71— l,n2 — 1 



(p) 



n — l,n„ — 1 



Therefore the type I MOPs uniquely exist if and only if 



(2.6) 



DfT :=detM,-i ^ 



(p) 



n—l.rLjy — l 



^0. 



The linear system arising from the type II conditions (|2.1[) has a matrix which 
is the transpose of (|2.5[) . Therefore the non- vanishing of the determinant (|2.6[) also 
guarantees the existence and uniqueness of the type II MOP. 

Suppose Dfi ^ 0. Then it is easy to see that the type II MOP has the deter- 
minantal formula 



(2.7) 



Pn{x) 



rr(l) rr(2} 

n.ni— 1 n. 712 — 1 



(p) 



Indeed, the right-hand side of (|2.7p is a monic polynomial of degree n. If we 
multiply the right-hand side of (|2.7p by x^Wj{x) and integrate with respect to x, 
we can perform these operations in the last column to obtain a determinant with 
two equal columns if fc < — 1. This proves the type II orthogonality conditions 

m- 

The type I MOPs have a similar determinantal expression. For j = 1, . . . ,p we 
have 



(2.8) J^{x)^—x 



1 



(1) 



n— 2,ni — 1 





H. 



0-1) 

n—2.7ij-i — 1 



H. 



n — 2 , n ^ — 1 



Ix-- - 



(j+i) 

n—2,nj-\-i —1 





(p) 

n — 2,np — 1 





These and similar determinantal formulas have recently been considered from 
the point of view of integrable systems in [H 111] . 
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2.3. Multiple integral representation. For what follows it is convenient to 
write 



and to introduce two sequences of functions /i , . . . , /^j and , . . . , g^j by 

(2.9) f,{x)^x^-\ 
and 

(2.10) gi+N,_dx) = x^'^Wjix), i = l,...,nj, j = l,...,p. 
Then the block Hankcl matrix (12.51) can be written as 

(2.11) 

and 

Dfi = det M.n = det [/^ fj{x)gk{x) dx]^^^^^^ 
For general m and n = \n\ we also write 

(2-12) ^W=K.^],.i....,™...i....,„, 
so that we have by (|2.7p 



Mfi = [mj,fc] „ , "ij.fe = / fj{x)gk{x) dx 



(2.13) 



and by jlH) and ((TTO)) 



Ps(a;) 



(2.14) Q.(a:) = ^4)(a:)t.,(x) = -i- 



51 (a;) ff2(a;) • • ■ gn{x) 



The following lemma is standard, see e.g. [45|, Proposition 2.10] where it is 
called a generalized Cauchy-Binet identity. 



Lemma 2.1. We have 

i*oo poo 



(2.15) Dn = -J ■ ■■ det [Mxk)]^^,^,_^ ■ det [g,{xk)]^,,^,_^ J] d^^' 
Proof. Expanding the two determinants on the right-hand side of (|2.15[) we 

get 

n 

det [/,(^^)],,, ■ det [gjixk]]^^, ^ ^^(_i)sg„.+sg„. -Q 

o- r k=l 

where the sums are for a and r over the symmetric group 5„. By (|2.1ip the 

right-hand side of (|2.15p is equal to 

(2.16) 



'™(TOr-i(fc),fc. 
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For any fixed a, we have that got ^ runs through 5„ as t runs through Sn- Hence 

n 

(2.17) ^(_i)sg„(.o.-) -Q ^ detMs = Dn. 

T k=l 

The equahty ||TT5| follows from and (PTZ)) . □ 

There is a similar multiple integral representation for the type II MOPs, which 
was stated for a special case in |14) , see also |32| . We emphasize that it is important 
here that fj{x) = x^~^ . 

Proposition 2.2. Assume 0- Then the type II MOP has the multiple 

integral representation 

(2.18) Ps(z) = _^x 

■■■ Y[{z~xk)-det [fj{xk)]^ „ • det [gj{xk)]^ „ J] '^^fc' 

Proof. Since fj{x) = x^^^ we have that det [fj{xk)] is a Vandermonde deter- 
minant, and therefore 

n 

Y[{z- Xk) ■ det [fji^k)]^ ,^,^ = det [fji^k)]^^,^,^^^^ ,^^, 

k=l 

where we have put 

fri+iix) = x", and Xn+l = z. 
Thus, by expanding the determinant we have 



k=i o-eS„+i fc=i 

and similarly 

n 

det [gA^k)]^^,^,_„ - J2 i-^r''^I[9rik)ixk). 

T&S„ k=l 

Integrating the product of the two above expressions with respect to Xi, . . . , a::„ 
we obtain 

(2.19) 

/OO poo Jf' 

■■■ \[{z-Xk)-Aet [/j(a;fe)]j- • det [gjixk)] ^j^^^^^^ \{ dxk 

n 

(TGSn + lTGSn k^l 

11 

TeS„ o-es„+i fe=i 

where we used the definition of rrij^k as given in p. lip also for j = n + 1. 

For each fixed t € Sn have that the sum over a in (|2.19p is equal to the 
determinant in the right-hand side of p.lSp and the proposition follows. □ 
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In an analogous way we find the following multiple integral representation for 
the linear form of type I MOPs, which is due to Desrosiers and Forrester 



Proposition 2.3. Assume Da ^ 0. Then the linear form of type I MOPs 
satisfies 



(2.20) 



dx — 



1 



k=l 



/oo poo ^ ^ 
-oo J —oo , 1 

Proof. Here we use the property 

n 

Y[{z- Xk)-^ ■ det [fj{xk)]^ ,^,^ ,^ = det [/,(a:^fe)]^- ,^i^...,„+i 



fe=i 

where now we put 



fn+lix) 



1 



and 



Xn+l = Z. 



Z — X 

The rest of the proof follows along the same lines as the proof of Proposition 12.21 
We omit the details, see also 1321. □ 



3. MOP ensembles 

3.1. Probabilistic interpretation. The multiple integral representations (|2.15p . 
p.lSp and (|2.20p have a natural probabilistic interpretation in case the product of 
determinants 

det [/^(a^fc)]^. fc^i ... „ • det ... „ 

is of a fixed sign for (xi, . . . , a;„) e R". That is, if it is always > or always < 0. 
Indeed, in that case it follows by (j2.15p that 

(3.1) V{xi,...,Xn) = ^det [/j(xfc)] . ■^et[gj{xk)].,,^i 

Z/.^ 115 111 

is a probability density function on R", where 

(3.2) Zn - Dr^nl 

is the normalizing constant (also called partition function in statistical mechanics 
literature), so that J ■ ■ ■ J 'Pixi, . . . , Xn)dxi ■ ■ ■ dxn — 1- 

The multiple integral representations (|2.18p and (|2.20p then show that 



(3.3) 



(3.4) 



Pn{z) = E 



,fc=l 



z e 



Yiiz-XkY 



.k=l 



z e C\] 



where the mathematical expectation is taken with respect to the p.d.f. (j3.ip . 

Thus Pn{z) is the average of the polynomials nfc=i('^ ^ ^k) where the roots 
distributed according to (|3.ip . In cases where the distribution p.ip 
can be interpreted as the eigenvalue distribution of a random matrix ensembles, 
one would call Pa the average characteristic polynomial. 
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3.2. Biorthogonal ensembles. A biorthogonal ensemble, see jl8j . is a prob- 
ability density function on K" of the form ()3.1|) with certain given functions /i, . . . , /„ 
and f/i, . . . ,g„, not necessarily of the form (|2.9p and (j2.10p . The p.d.f. is invariant 
under permutations of variables. We think of the ensemble as giving us n random 
points or particles Xj on the real line, and so it is a random point process. 

A biorthogonal ensemble is a special case of a determinantal point process, see 
e.g. [45^ [63) . This means that there is a correlation kernel Kn{x, y) so that 



1 



and so that marginal densities (m point correlation functions) are determinants 

{n — m)\ 



/OO pC 



Vixi, . . . ,Xn)dx> 



m+1 



' ■ dXr, 



■ det [Kn{xj,Xk)] 



l,...,m 



n—rn times 



Taking for example m — 1 we have that ^Kn{x,x) is the mean density of points, 
that is 



1 r 

— / Kn{x,x)dx 

n J a 



is the expected fraction of points lying in the interval [a, h]. 

In a biorthogonal ensemble, the correlation kernel can be written as a bordered 
determinant 



(3.5) 



-1 



det M„ 



M„ : 

fn{x) 

9i{y) ■ ■ ■ 9n{y) 



where M„ is the matrix 

M„ = [mj- fc] . ^^^^ 



fj (x)gk {x) dx 



In the formulation of the biorthogonal ensemble p.ip . we have some freedom 
in choosing the functions /i,...,/„. and gi, . . . , gn- Indeed, if (pi, . . . , 4>n and 
rpi, . . . jipn are functions with the same linear span as the /j's and ff^'s, respectively, 
then we could use these functions instead. A particular nice form appears if the 
functions and ^pk are biorthogonal, i.e., 



Then the representation 



(3.6) 



(j)j{x)ilJk{x)dx = Sj^k 
reduces to 

(j>i{x) 
In : 

V'i(y) ■ • • i'niy) 
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3.3. OP ensembles. If fj{x) ~ gj{x) — x^ ^ ^Jw{x), j — 1, . . . ,n for some 
non-negative weight function w on M, then 

1 1 " 

(3-7) — det [fj{xk)] ■ det [gj{xk)] = — IT ^^'^ ~ ^^^^ ' 11 ^^^''^ 

" " l<j<fc<ri fc=l 

which is indeed of constant sign. This is the form of the joint p.d.f. for the eigen- 
values of a unitary random matrix ensemble 

■^cxp(-TTV{H))dH, w(x)=e-^(^', 

defined on n x n Hermitian matrices H, see j29j . 

In this case the biorthogonal functions take the form 

(3.8) 0j(x) ='0j(x) =p-j^i{x)\/w{x) 

where Pj-i is the orthonormal polynomial of degree j — I with respect to the weight 
w on M, and by (|3.6p , 

Kn{x,y) = ^ w{x) \/w{y)^ Pj {x)pj (y) 

3=0 

is the correlation kernel, which in this situation is also called the OP kernel or the 
Christoffel-Darboux kernel. 

The OPs are characterized by a 2 x 2 matrix valued Riemann-Hilbert problem 
due to Fokas, Its, and Kitaev [39) . 

• y : C \ R ^ C2^2 is analytic, 

. Y+{x)^Y^{x)(l ""f^y xeR, 

• y(z) = (/2 + C'(l/z))diag(z" asz->oo. 

The correlation kernel for the OP ensemble can be given directly in terms of the 
solution Y of the RH problem 



(3.9) Kr,{x,y)^^-^-—^^/^^/^){0 l)Y+iy)-'Y+ix) 

This follows from an explicit formula for Y in terms of the orthogonal polynomials 
Pn and Pn-i, and the Christoffel-Darboux formula for orthogonal polynomials. 

3.4. MOP ensembles. We have a MOP ensemble if /i, . . . , /„ and gi, . . . ,gn 
are given by (I2J1) and (fTTO]) . and if 

(3.10) det [Mxk)]^^,^,_^ ■ det [gA^k)]^^,^,_^ 

has constant sign. The case p = 1 reduces essentially to the OP case. 

For a MOP ensemble we have that the correlation kernel Kn given by the 
determinant ()3.5p has another expression in terms of the RH problem for multiple 
orthogonal polynomials. MOPs (with p weights) satisfy a (p -I- 1) x (p + 1) matrix 
valued RH problem 68 

• y : C \ R ^ c{p+-L)>^{p+i) is analytic. 
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(\ Wi W2 

1 

■•• 



Wp\ 





X e 



■■• 

Vo 1/ 

• Y{z) = (/p+i + 0(l/z)) diag (z" z""! • • • z""") as z oo. 

The correlation kernel for the MOP ensemble is given as follows in terms of the 
solution Y of the RH problem 

/1\ 



(3.11) K^{x,y)^ 



1 



2'Ki{x — y) 



(0 w,{y) ■■■ Wp{y))Y+{y)-'Y+{x) 







The proof is based on a Christoffel-Darboux formula for MOPs, see |14| for the 
case p = 2 and 25j for general p. An extension to MOP of mixed type is given in 



4. Special classes 



The condition that p.lOp is of constant sign is automatically satisfied in the 
OP case, but it becomes relevant in the MOP case. We call it the MOP ensemble 
condition. It is of interest to identify classes for which the MOP ensemble condition 
holds. In the literature on Hermite-Pade approximation a number of special classes 
of MOPs were identified for which rather complete convergence results could be 
established. These are in particular the Angelesco systems [2], 140] and the Nikishin 
systems, see e.g. [H 1201 1341 142L 150] . It turns out that for these special classes the 
MOP ensemble condition holds. Before we turn to that, we make some preliminary 
observations. 

4.1. Preliminary observations. The first observation is that the product of 
determinants (|3.10[) is invariant under permutations of the x^'s. It is also clear that 
p.lOp is zero in case two or more of the Xk 's coincide. Therefore we may restrict 
ourselves to strictly ordered sets of points 

(4.1) Xi < X2 < ■ ■ ■ < Xn- 

The second observation is that the first factor in (|3.10p is a Vandermonde determi- 
nant (due to the fact that fj{x) = x^^^) 

det [fjixk)]- ,^^^ ^ ,^ = Ylixk- Xj) 

■j<k 

which is positive for ordered points (14. ip . Therefore the MOP ensemble condition 
comes down to the condition stated in the following lemma. 

Lemma 4.1. The MOP ensemble condition is satisfied if and only if either 

(4.2) det[<?j(xfe)]^.^^^_ „>0 
whenever xi < X2 < ■ ■ ■ < Xn, or 

(4.3) det[g,(xfe)]^.^^^ <0 
whenever xi < X2 < ■ ■ ■ < Xn- 
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Since Vandermonde-like products will appear frequently in what follows we use 
the abbreviations 

(4.4) A{X)= H (xk-x,) 

l<j<k<n 

if X = (xi, . . . , Xn), and 

n m 

(4.5) A{X,Y)^l[Y[{xk-yj) 

k=ij=i 

if X = (xi, ...,Xn) and Y = (yi, . . .,ym)- 



4.2. Angelesco ensemble. The weights wi, . . . ,Wp are an Angelesco system 
if there are disjoint intervals Fi, . . . , Fp C M, such that 

supp(wj)cFj, j^l,...,p. 

We write F^ — [aj, Pj] and without loss of generality we assume that 

(4.6) Pj < aj+i , for j = 1, . . . , p - 1. 

We extend Wj to all of R by defining wj (x) = if a; G M \ Fj . 

An Angelesco system always gives rise to a MOP ensemble. Indeed it is easy 
to see that in the Angelesco case det [^^(xfe)] is of block form, and it can only be 
non-zero if rij of the points Xk belong to F-,- for j — I, . . . ,p, and so this is what we 
will assume. Then writing Nj = = 0, we have that 



XNj^i+k & Fj 



1,.. 



Because of the orderings (|4.ip and (|4.6p the determinant then has a block diagonal 
form where the ith block is 



/ w,{x^;^) 

Xl 'Wi[Xl') 



Wi{Xn}) 



V(4")' 

whose determinant is 

n ( 



Wi [x\ ') 



Xl — x)' 



l<j<k<ni k=l 



k=l 



where A^*' = {x^^' , . . . , x^nj). The resuh is that 
(4.7) det [g,{xk)] = H ^(^''^) ' 11 ^^(4^^) 



fc=i 



and this is > for every choice of xi < • ■ • < x„. 

Thus, by Lemma (|4.ip . an Angelesco system gives rise to a MOP ensemble and 
we call it an Angelesco ensemble. We see from the above calculation that the joint 
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p.d.f. in an Angelesco ensemble is 



1 

p n 



(4.8) —det[fj{xk)]det[gj{xk 



" 1=1 l<i<j<p 1=1 fe=l 

4.3. AT ensemble. Assume that wi,. . . are weights defined on a fixed 
interval F C M. Then wi, . . . ,Wp are an AT system on F if the functions gj are 
an algebraic Chebyshev system on F. This means that every non-trivial linear 
combination 

n 

has at most n — 1 zeros in F. Equivalently, an algebraic Chebyshev system means 
that 

(4.9) det [g,{xk)\ + 

for every choice of distinct points in F. The property of being an AT system 
also depends on the multi-indices ni, . . . , rip. 

If the weights Wj in an AT system are continuous functions on F, then it clearly 
follows from (|4.9p by continuity that det [^^(a;*:)] has constant sign (either > or 
< 0) whenever the are strictly ordered points in F. Therefore, in that case, we 
have a MOP ensemble by Lemma l4Tl which we we will call an AT ensemble. 

4.4. Nikishin ensemble. 

4.4.1. Definition of a Nikishin system. Certain AT systems with special prop- 
erties were first described by Nikishin [58| and are therefore called Nikishin systems. 
We state it first for p = 2 continuous weight functions wi , W2 defined on an interval 
Fi c R. 

The Nikishin assumption is that the ratio 1x12/1111 can be written as a Markov 
function for a non-negative weight function (or more generally a measure) supported 
on an interval F2, disjoint from Fi, that is, if 

(4.10) ^^^j 

where v{s) is a non- negative weight function with 

supp(i;)=F2, F2nFi=0. 

We choose the + sign in (|4.10|) if F2 lies to the left of Fi; otherwise we choose the 
— sign. Then we call wi,W2 a Nikishin system on Fi for the intervals Fi,F2. 

A Nikishin system withp > 3 weights is defined inductively. Suppose supp(wj) — 
Fi for all J — 1, ... ,p, where Fi is an interval. Suppose 



(4.11) -^^± ^^dy, xeFi, j^2,...,p, 

where F2 n Fi =0 and where V2, ■ ■ ■ ,Vp is a Nikishin system on F2 for the intervals 
F2,...,Fp. Then we call wi,...,Wp a Nikishin system on Fi for the intervals 
Fi, . . . ,Fp. 
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Note that in a Nikishin system two consecutive intervals Tj and Fj+i are dis- 
joint. However, if |j — fc| > 2, then Tj and Tk may very well have a non-empty 
intersection. 

This construction might not seem very natural at first sight, but it is actually 
a very beautiful structure. A main result is that for multi-indices n = (ni, . . . , Up) 
such that 

Tij > rij+i - 1, j = 1, . . . ,p- 1 

a Nikishin system is an AT system, see |59j . and therefore the type I and type II 
MOPs exist. As we have seen in the previous subsection, there is also an associated 
MOP ensemble, which we call a Nikishin ensemble. 

4.4.2. Nikishin ensemble with 2 weights. Here we show that a Nikishin ensemble 
has a natural interpretation as the marginal distribution of an extended ensemble. 
The following calculations are due to Coussement and Van Assche |22j . 

For reasons of clarity we take p = 2 and we assume that r2 is to the left of Fi. 



Then for xi, 



, Xn in Fi we have 
wi{xi) 

XiWi{xi) 



det [gjixk)] = 



W2ixi) 



'^W2ixi) 



Y[ ^"i(a;fe) 

k=l 



Wi{x2) 



W2{X2) 



1 

Xi 



ni — 1 

^1 



K2-1 W2(xi) 

1 wi{xi) 



Wi{Xn) 
XnWi{x„) 

W2ix„) 



W2{X2) 
Wi(x2) 



W2{x„) 

■i«i(a;„) 



712-1 W2{X„ 

" wi (x„ 



Now we replace each ratio by the integral (|4.10p , we use t/j as the integration 

variable in row ni -I- j, and we take the integrals as well as the factors v{yj) out of 
the determinant, to obtain 



n wi{xk) 



k=l 




Since 



ni > n2 — 1 



1 


1 


1 


ni — 1 

^1 

1 


1 


1 


xi-yi 
xi 


x2-yi 


Xn-yi 

Xn 


xi-y2 




Xn-y2 


"2-1 

^1 






xi-y„2 




x„-yn2 
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we can perforin elementary row operations to reduce the remaining determinant to 



1 


1 


1 




1 


1 


1 


ni — 1 

1 


1 


„Tli-l 

■''n 

1 


"2 

=n^r 


rii —1 

^1 

1 


1 


T,rii-1 
1 




x2-yi 


x„-yi 
y2 


xi-yi 


2^2-yi 


a;7i-yi 


X1-V2 
"2-1 




Xn-y2 

"2-1 

y2 




xi-y2 
1 




x„-y2 
1 


Xl-yn2 




X„-yn2 




Xl-yn2 




x„-y„2 



which is a mixture of a Vandermonde and a Cauchy determinant. It can be evalu- 
ated to give 

4^ A{X)A{Y) 

where X = (xi, . . . ,a;„) and F = (j/i, . . . jj/nj. Thus 
(4.12) det[gjixk)] 

n n n "2 "2 A('"V^A 



Now we symmetrize the multiple integral with respect to the integration variables 
Uj (which is a standard trick in determinantal point processes). That is, for any 
permutation a G Sn2 we make the change of variables yj i-^ Uaij) and we average 
over all permutations a in Sn2 ■ Using the fact that 

"2 

E (-1)""'^ n y^' = [yr^],..,...,„, = A(n 

we then obtain that (I4.12p is equal to 
(4.13) det[gjixk)] 



1 " f f A(Y]'^ 

- n («.) A(x, . . ^ /^^ n "(y,)^, n 

A; — 1 J — -L J — -L 



■'^A(x,y)-ii"^^- 



The joint p.d.f. for the Nikishin ensemble is therefore (since det [fj{xk) 
A{X)) 



(4.14) — det[/,(xfe)]det[5,(a:fe)] 
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By dropping the integrals over the yj variables, we can view (|4.14p as a marginal 
density of an extended ensemble defined by the joint p.d.f. 

(4.15) Vext{xi,. . .,Xn,yi,- ■ ■ ,2/712) 

= z::;^n n -(2/.) • A(x,y) 

defined for xi, . . . , a;„ G Fi and yi, . . . , yn2 G ^2- Note that the factor A{X, Y) in 
(|4.15p is positive, since r2 lies to the left of Fi so that Xk > yj for every k — 1, . . . , n 
and j = 1 , . . . , ?i2 . 

4.4.3. Nikishin ensemble with p > 2 weights. The above considerations can be 
extended to general p > 2. Let wi, . . . , Wp be a Nikishin system with p weights for 
the intervals Fi, . . . , Fp. Assume that 

rij > Tij+i - 1, for j = 1, . . . ,p - 1. 

Then the joint p.d.f. for the Nikishin ensemble is a marginal density of an extended 
ensemble defined on 

p 

Ff 1 X F^^ X • • • X F^p, where = ^ n„ 
with joint p.d.f. of the form 

(4.16) i fr TT ^"-'^ i^) ■ n^-i^(^^'^)' 

^^/JH ^ ' ^ n^zjA(xa),xa+i)) 

where w^^^ is a certain weight function on F^ for j = 1, . . . ,p, with iv^^^ ~ wi. Here 

X« = 4] )eFf 

and Zn is a normalizing constant. 

5. Weak asymptotics 

An important question about a sequence of polynomials with increasing degrees, 
is about the asymptotic behavior as the degree tends to 00. 

5.1. Vector equilibrium problems. To describe the weak asymptotics of 
the type II MOPs, as well as the convergence for the Hermite-Pade rational approx- 
imation problems, vector equilibrium problems were identified that are relevant for 
the Angelesco and Nikishin systems. Here we show that these equilibrium prob- 
lems have a natural interpretation in terms of the joint p.d.f. 's of the Angelesco and 
Nikishin ensembles. 

We assume that we are considering MOPs Ph for a sequence of multi-indices n 
such that n = \n\ 00 and Uj — > 00 for every j = 1, . . . ,p in such a way that 

TL ■ 

(5.1) —^rj for j = l,...,p. 

n 

The limiting ratios Vj should satisfy 

p 

0<rj<l, ^r,. 
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Here and in the following we use 

I{f^)= //log , ^ , dfi{x)dfi{y) 
J J \x-y\ 

to denote the logarithmic energy of the measure ^, and 

lif^,!^) =11 log] — ^ — -dn{x)dv{y), 
J J \x-y\ 

which is the mutual logarithmic energy of the two measures /i and v. 
For a discrete measure we introduce the reduced logartithmic energy 

I*{^J^) = jj ^og j-^—^dn{x)dn{y) 

Note that, if 

n 
k=l 

is the point counting measure of X = (xi, . . . , x„) E R", then 

\ogA{Xf = -r{i,x). 

5.2. Angelesco system. The Angelesco ensemble has the joint p.d.f. 
A configuration of points X in an Angelesco ensemble is of the form 

(5.2) A = (A«,...,A(p)), where A« = . . . , x.^) £ 
The most likely configuration minimizes 

(5.3) -^logA(A0))2-5] logA(AW,A(^)) + ^^Q,(x«) 

j=i i=i j=i k=i 

where Wj — e"*^^ , among all X of the form (j5.2p . Introducing the normalized point 
counting measures 

(5.4) Vj^-vxu), j = l,...,p, 

n 

we can rewrite ()5.3p . after dividing by n^, as 

(5.5) E^^C'^.) +E E ^(^-^^) + ^E/ Q^i^)d^A^)- 

j=l 1=1 j=i+l j=l 

In the limit (|5.ip we forget about the discreteness of the measures I'j. Then 
instead of minimizing (|5.5p among all vectors of measures (z^i, . . . , Vp) with Vj a 
measure of total mass rij/n on Tj of the form (j5.4p . we come to minimize the 
energy functional 

p p~i p 

(5.6) £:(/!!, Ai2, Mp) =E^('-'j)+E E ^(t^O'f^k) 

J=l J=l fc=J+l 

among all vectors of measures {fii, . . . , /Zp) with 

(5.7) supp(/ij) C Tj and y d/ij = rj. 

Under the assumption that each Tj is compact, and that Wj{x) > almost 
everywhere on Tj, Gonchar and Rakhmanov [40j showed that the zeros of the 
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type II MOP Pfi are distributed according to the minimizer (/ii, //2, • • • , Mr) for 
the vector equihbrium problem. More precisely, for every j — 1, . . . there are 
rij simple zeros of in Tj, say x'f\ . . . ,Xnj, and the normalized zero counting 
measure 



n. ^ ^ 



n 

k=l 

converges in the limit (|5.1|) weakly to fij. 

Under the same conditions, it seems likely that the vector of normalized count- 
ing measures {vi, . . . , Vp) (as in (j5.4p ') of a random point (j5.2p from an Angelesco 
ensemble tends to the vector of nonrandom measures (^i, . ■ . ,Mp)i almost surely, 
but this has not been established rigorously. 

In a situation of varying weights in an Angelesco ensemble, such as for example 

«;,(a;)=e-"^^(-), x&T,, j = l,...,p, 

we have to add external field terms to (|5.6p and the relevant energy functional 
becomes 

p p-i p ^ r 

(5.8) E^(^:') + E E n^^J,^^k) + J2 / vA^)d^i,{x). 

3 = 1 3 = 1 k=3 + l 3 = 1 

This concept is well known in the orthogonal polynomial case, see |62) for the 
standard reference on logarithmic potential theory with external fields. 

5.3. Nikishin system. Similar considerations apply to the joint p.d.f. (|4.16p 

in the extended Nikishin ensemble. Here the most likely configuration 

where 

X(^)^(4^\4^\...,a;(^])Grf 
is the one that minimizes 

p p—l p Nj 

(5.9) - E A(X(^) )2 + J2 log A(X(^) , ^(■'■+1) ) + E E (^fe'^ ) 

j=i 3=1 3=1 k=i 

where w^^^ = e"*^^. In terms of the normalized point counting measures 

(5.10) = j = l,...,p, 

n 

the expression (|5.9p is, after dividing by n^, 

p p — 1 p 

(5.11) E^*('^^) -E^(^^'^^+i) + ^E / QA^)d'^ji^)- 

j=i j=i ^ i=i 

The measure Vj has total mass 



/ ai^3 ^ — = - > n^, 
J n n^. 



and in the limit (|5.ip we have 



Tt ^ ^ 
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So in the limit (|5.1|) we are led to the following energy functional 



p-i 



i=i i=i 



(5.12) 

among all vectors of measures {fii , . . . , 
(5.13) supp(^j) C Tj 



lip) with 

and / dfij — i 



In the special case that rii = n2 
and then the normalizations are 



d/ij — 1 



rip we have that all are equal to 1/p, 



P 



,P- 



5.4. Interaction matrix. In both an Angelesco system and a Nikishin system 
we minimize 

p V 
3 = 1 k=l 

with a certain positive definite interaction matrix C — (cjk)- 
The Angelesco interaction matrix is a full matrix 

1 ... 1 i\ 

? ? V 

2 2 2 



/I 

1 
2 



(5.14) 



1 1 

2 2 



1 



1 1 

! f 

2 2 



2 

1/ 



and the Nikishin interaction matrix is a tridiagonal matrix 



/ 1 



(5.15) 



C 







V 









_1 

2 

1 / 



The Angelesco interaction matrix (|5.14p shows the repulsion that exists between 
fj,j and ^fc when j ^ k. However, the repulsion is only half as strong as the inner 
repulsion between each measure fij itself. 

Since the non-zero off-diagonal entries in (j5.15p are negative, there is an at- 
traction in a Nikishin ensemble between two consecutive measures fij and Hj+i for 
j = 1, ... .J) — 1. The measures fij and fik with \j — fc| > 2 do not interact. The 
Nikishin interaction also arises in the asymptotic analysis of eigenvalues of banded 
Toeplitz matrices |35j as well as in the two-matrix model [36j where it appears 
with both an external field and an upper constraint. 
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